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Abstract. We introduce a simplified model for the minimization of the elastic energy in thin shells. This 
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> 
O 



o 



model is not obtained by an asymptotic analysis. The nonlinear simplified model admits always minimizers 
by contrast with the original one. We show the relevance of our approach by proving that the rescaled 



minimum of the simplified model and the rescaled infimum of the full model have the same limit as the 



thickness tends to 0. The simplified energy can be expressed as a functional acting over fields defined on the 



mid-surface of the shell and where the thickness remains as a parameter. 
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1. Introduction 



This paper is devoted to introduce and justify a simplified model for nonlinear clastic shells. Let w be a 
£\j , bounded Lipschitz domain of M 2 and be a smooth function from oJ into K 3 (see the detailed assumptions 

on <p in Section 2) and set S — 4>{oS). We denote by n an unit vector field normal to S and by <fr the 
\ map (si,S2,S3) — > (f>(si,S2) + S3n(si,S2)- The elastic shell is defined by Qs — $(wx] — S,8[) and we 

consider that it is clamped on a part of its lateral boundary Tq,s = $(70 x] — 5, 5[), where 70 C du) . The 
energy density is denoted W and we assume that Qs is submitted to applied body forces f K s whose order 
with respect to S depends upon a parameter k (see the order of f K .s below). The total energy is given by 
J KtS (v) = [ W{E[v)) - I f K . s ■ (v - I d ) if det(Vw) > and where E(v) = l/2((Vi;) T Vu - I3) is the 
Green-St Venant's tensor and Id is the identity map. We set 

• 1— 1 

X 

5h ' m K ,s = inf J K ,s(v), 

where Us is the set of admissible deformations (which are equal to the identity map on Tq^s). The Korn's 
type inequalities established in [6] (see also [f 2]) allow us to prove that if the order of f Kl s is equal to S 2k ~ 2 
for f < k < 2 (or S K for k > 2), then the order of m K ,s is i5 2k_1 . 

Even for a classical St-Venant-Kirchhoff's material, proving the existence of a minimizer for J Kt s is still 
an open problem. The aim of this paper is to replace the above minimization problem by a minimization 
problem for a simplified functional s defined on a new set B><5,7 and which admits a minimum 

v£bj i10 

of the same order as m Kt s- This approximation is justified if one shows that 

m K ,s - m% s 
hm — — : — - = 0. 
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In the present paper we show this result in the case n = 2 (other critical cases will be investigated in 
forthcoming papers). 

The expression of J* $ and the choice of B,5 j7o rely on the decomposition technique introduced in [6] . 
Let us recall that a deformation v of the shell Qs, whose "geometrical energy" -||dist(Vv, SO(3))||x,2(q 4 )- is 
at most of order S 3 ^ 2 , is decomposed as (see [6] or Theorem 3.1 below) 

v{x) = V(si, s 2 ) + s 3 R(si, s 2 )n(si, s 2 ) + v(si,s 2 ,s 3 ), x = $(s), for a.e. s = (s 1 ,s 2 , s 3 ) E ujx] - 5, 6[. 

The field V stands for the mid-surface deformation, the matrix field R takes its values in SO(3) and represents 
the rotations of the fibers and v is the warping of theses fibers. It is also shown in [6] that the fields V, 
R and v satisfy the natural boundary conditions on 70 and on 70 x] — S, S[ and that they are estimated in 
terms of \\dist(Vv, SO(3))\\i J 2fQ 5 ^ and 6. With the help of these estimates, we justify the simplification of 

the Green-St Venant's strain tensor E(v) in order to give a simplified matrix E(v) which depends on the 

_ . . dv 

triplet v = (V, R, v) associated to a deformation v. This matrix depends linearly upon — — and on the first 

os 3 _ 

partial derivatives of V and R (see Section 5) but which is nonlinear with respect to (V, R, v). 

Then we define the set 1D>5, 70 of admissible triplets v = (V, R, v) and we derive the simplified total energy 

J* j(5 (v) as follows. Firstly we replace f Q W(E(vj) by Q(E(v)) where Q is a quadratic form which is 

assumed to approximate W near the origin. Secondly we add two penalization terms in order to approach 

dV 

the usual limit kinematic condition — — = Rt Q and to insure the coerciveness of J% Finally in the term 

ds a 

involving the forces we neglect the contribution of the warping v. As announced above we prove that J* s 
admit minimizers on B,5 j7o . We justify the approximation process described above in the case k = 2. 

In some sense, the introduction of J* s can be seen as a nonlinear version of the approach which leads to 
the simplified Timoshcnko's model for rods, the Reisner-Mindlin's model for plates and the Koiter's model 
for shells in linear elasticity. 

As general references on the theory of nonlinear elasticity, we refer to [8] and [24] and to the extensive 
bibliographies of these works. A general theory for the existence of minimizers of nonlinear elastic energies 
can be found in [1]. For the justification of plate or shell models in nonlinear elasticity we refer to [9], [10], 
[11], [13], [15], [18], [23], [25], [26]. The derivation of limit energies for thin domains using T-convergence 
arguments are developed in [14], [15], [22], [23]. The decomposition of the deformations in thin structures is 
introduced in [17], [18] and a few applications to the junctions of multi-structures and homogenization are 
given in [2], [3], [4]. The justification of simplified models for rods and plates in linear elasticity, based on 
a decomposition technique of the displacement, is presented in [19], [20]. In this linear case, error estimates 
between the solution of the initial model and the one of the simplified model are also established. In some 
sense, these works give a mathematical justification of Timoshenko's model for rods and Reisner-Mindlin's 
model for plates. 

The paper is organized as follows. Section 2 is devoted to describe the geometry of the shell and to 
give a few notations. In Section 3 we recall the results of [6]: decomposition of a deformation of a thin 
shell, estimates on the terms of this decomposition and two nonlinear Korn's type inequalities. Section 4 is 
concerned with a standard rescaling. We present the simplification of the Green-St Venant's strain tensor of 
a deformation in Section 5. We also introduce the set B>i, 7o of admissible triplets v = (V, R, V) and we prove 
Korn's type inequalities for the elements of B,5 j7o (see Corollary 5.3). In Section 6 we consider nonlinear 
clastic shells and we use the results of [6] to scale the applied forces in order to obtain a priori estimates on 
m K fi. Section 7 is devoted to introduce the simplified energy J S K s and to prove the existence of minimizers. 
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In Sections 8 and 9, we restrict the analysis to k = 2. We prove that 



,. m 2 ,s .. ™> 2 s 

lim — -r— = hm — -— = m 2 

6-K) S 3 <5^o S 3 

where mf is the minimum of a functional defined over a set of triplets. In Section 10, we give an alternative 
formulation of the minimization problem for J* s through elimination of the variable v. Then we obtain 
that m 3 K s is the minimum of a functional which depends only upon (V,R). At last an appendix contains 
an approximation result for the elements of I\ 7o and an algebraic elimination process for quadratic forms. 
The results of this paper were announced in [7] . 

2. The geometry and notations. 

Let us introduce a few notations and definitions concerning the geometry of the shell. 

Let w be a bounded domain in R 2 with lipschitzian boundary and let <j) be an injective mapping from 

deb 

to into R 3 of class C 2 . We denote S the surface 4>{Uj). We assume that the two vectors — — (si,s 2 ) and 

OSl 

o — i s ii s 2) are linearly independent at each point (s\, s 2 ) G oJ. 
We set 

osi os 2 ||tiAt 2 || 2 

The vectors ti and t 2 are tangential vectors to the surface S and the vector n is a unit normal vector to this 
surface. We set 

Clg = ljx] - S, S[. 
Now we consider the mapping $ : w x E — > M 3 defined by 

(2.2) $ : (s 1 ,s 2 ,s 3 )i — > x = (j}(s 1 ,s 2 ) + s 3 n(s 1 ,s 2 ). 

There exists d €E (0, 1] depending only on S, such that the restriction of $ to the compact set $}g = 
ZJ x [— 60, 5q\ is a C 1 — diffcomorphism of that set onto its range (see e.g. [21]). Hence, there exist two 
constants cq > and c\ > Co, which depend only on <p, such that 

(2.3) Vsen 5 „, c < |||V»$(a))||| < c u and for x = $(a) c < || | V«* _1 (a;))|| j < c x . 
Definition 2.1. For 5 e (0,(5 ], the shell Qg is defined as follows: 

Qs - $(n 4 ). 

The mid-surface of the shell is S. The fibers of the shell are the segments <&({(si, s 2 )}x] — S, S[), (si,s 2 ) G u>. 
The lateral boundary of the shell is Tg = $(<9a;x] — 5,8[). In the following sections the shell will be fixed on 
a part of its lateral boundary. Let 70 be an open subset of dui which made of a finite number of connected 
components (whose closure are disjoint). We assume that the shell is clamped on 

r 0l4 = $(7ox]-MD- 
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The admissible deformations v of the shell must then satisfy 
(2.4) v = I d on To.8 

where Id is the identity map of M . 

Notation. From now on we denote by c and C two positive generic constants which do not depend on 
5. We respectively note by x and s the generic points of Qg and of Q5. A field v defined on Qg can be also 
considered as a field defined on Clg that, as a convention, we will also denote by v. As far as the gradients 
of field v, say in (W 1 ' 1 (Q$)) 3 , are concerned we have V x v and \7 s v = \7 x v.S7& for a.e. x = <£>(s) and (2.3) 
shows that 

c|||V x ^)|||<|||V s ^)lll<C|l|V^)|||. 

3. Korn's type inequalities for shells. Decomposition of a deformation. 

We first recall the Korn's type inequalities for shells established in Section 4 of [6] . Let v be an admis- 
sible deformation belonging to (-ff 1 (Q,s)) 3 and satisfying the boundary condition (2.4). Setting V(si,s 2 ) = 
1 [ s 

— I v(s\,S2,t)dt a.e. (s\, s 2 ) € co, we have 
2o J-s 

j \\v- Id\\(L*(Q s ))3 + W^xv - h\\(Li(Q 5 )r < C(5 1/2 + \\dist(V x v, SO(3))\\ L 2 (Qs) ) , 
1 ||(« - I d ) - (V - 0)||( ia ( e4 ))s < C5(5 1 ' 2 + \\dist(V x v, S0(3))\\ LHQs) ), 



(3.1) 
and 
(3.2) 



C, 



v ~ Id\\(L*(Q 5 ))* + ||V X «- 1311(1,2(0)9 < j\\dist(V x v,SO(3))\\ L 2 (Qs) , 

(v - i d ) - (v - 4>)\\(L* { Q 6) r < c||dist(v^, so(3))\\ LHQs) . 

Inequalities (3.1) are better than those (3.2) if the order of the geometric energy \\dist(V x v, SO(3))\\l2(q s - ) 
is greater than S 3 ^ 2 . 

Now the theorem of decomposition of the deformations established in [6] (see Theorem 3.4 of Section 
3) is given below. 

Theorem 3.1. There exists a constant C(S) which depends only on the mid-surface of the shell such that 
for all deformation v belonging to (H 1 (Q$)) 3 and satisfying 



(3.3) 



\\dist(V x v,SO(3))\\ LHQs) <C(S)6 3 ' 



then, there exist V E (_ff 1 (w)) 3 , R e (iJ 1 (o;)) 3><3 satisfying H(si,s 2 ) € SO(3) for a.e. (s\,s 2 ) € u and v 
belonging to (ff 1 (Q < 5)) 3 such that for a.e. seDj 

(3.4) v(s) = V(si,s 2 ) + s 3 R(si,s 2 )n(si,s 2 ) +v(s), 

1 f s 

where we can choose V(si,s 2 ) = — r / v(si,s 2 ,t)dt a.e. (si,s 2 ) € uj, and such that the following estimates 

2<3 J-s 

hold: 

' INI(^ ( n 5 ))3 < C5\\dist(V x v,SO(3))\\ LHQs) 
||V s tJ|| (L 2 (n5 )) 9 < C\\dist{V x v,SO{3))\\ L 2 (Q5) 



(3.5) 



<9R 



ds a 
dV 
ds a 



< 



C 



(£ 2 M) 
Rt 



— \\dist{V x v,SO{3))\\ L i (Q5) 
C 

(L 2 (w)) 3 ~ (5 1 / 2 

(i 2 (^)) 9 



< — \\dist(V x v,SO(3))\\ L 2 iQs) 
< C\\dist(V x v,SO(3))\\ L 2 {Q5 y 



□ 



f 5 - 

Due to (3.4) and to the definition of V, the field v satisfies / v{s\, s 2 ,t)dt — a.e. (s\, s 2 ) G u. 

J-s 

If the deformation v as in Theorem 3.1 satisfies the boundary condition (2.4) then indeed 



(3.6) 



V 



on 70 . 



Moreover due to Lemma 4.1 of [6], we can choose R and v in Theorem 3.1 above such that 



(3.7) 



R = I 3 on 70, v = onr ,«5. 



From estimates (3.5) we also derive the following ones: 

C 



(3.8) 



\R--h\\{L*{ U ))» < lpj2 

c 



\\dist{V x v,SO{S))\\ 



L 2 (Qs) 



[ \\V - cj>\\ (LHui))3 < _||^(V^,50(3))|| L2(&) . 



4. Rescaling Q s 

As usual, we rescale £1$ using the operator 

(n s w)(s 1 ,s 2 ,S 3 ) = w(s 1 ,s 2 ,SS 3 ) for any (51,52,5*3) € Cl 

defined for e.g. w £ L 2 (Qs) for which (Ugw) e L 2 (Q). Let v be a deformation decomposed as (3.4), by 
transforming by Us we obtain 

U s (v)(si,S2,S 3 ) = V{s 1} s 2 ) + 5S 3 K(s 1 ,S2)n(s 1 ,s 2 ) + Us(v)(si,s 2 ,S 3 ), for a.e. (si, 52,5*3) G O. 

The estimates (3.5) of v transposed over f2 are (notice that Us ( ^ — ) = -r — — — ) 

\os 3 J oS 3 



( ||n 4 (t;)|| (i2(n)) . < C5 1 / 2 ||d^(V :c i;,50(3))|| L2(a 



(4.1) 



dU s (v) 



dsi 
dUs(v) 



C 



ds 2 
dUsiy) 



dS 3 



(i2(n))3 < ^\\dist(W x v,SO(3))\\ L2{Q5) 
( ^(Q))3 " ^IM^(V^,50(3))|U 2(&) 
^C^H^v^,^))!!^^), 

(i 2 (0)) J 



5. Simplification in the Green-St Venant's strain tensor. 

In this section we introduce a simplification of the Green-St Venant's strain tensor E(v) = 1 /2((Vw) T Vw— 13) . 
Let v be a deformation of the shell belonging to (H 1 (Qs)) 3 and satisfying the condition (3.3). We decompose 
v as (3.4). We have the identity 



(V x v) T V x v - I 3 = (V x v - R) T R + K T (V x v - R) + (V x w - R) T (V X « - R). 
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In order to compare the orders (of the norms) of the different terms in the above equality, we work in the 
fix domain f2 using the operator Us. Thanks to estimates (3.5) we get 

' |M(V,, - Rf R + R T (V x v - R))|| (Ll(n)) .x, < c \\^ xV ,SO m \ LHQe) 



(5-1) 



||Oj((V^ - R) T (V x v - R))|| (i i (n))Sx . < c 



s 1 / 2 

\\dist(V x v,SO(3))\\ L 2 {at y 2 



S 1 / 2 



In view of (3.1), these estimates show that the term Yls({V x v — R) T (S7 x v — R)) can be neglected in E{v). 
Now we have 



Then 



dv _ / dn , dv _ / dn 

^—=V x v(ti+s 3 — ), — = V x v (t 2 + S3-^— 

os\ asi' ds 2 0S2 



Il s{ V x v - R) (t a + 5S 3 ^-J = - Rt a ) + 5S 3 ^-n + ^ 



dv _ 
— = V x vn. 

ds 3 



dn 



dn dn 

First, we can neglect the term SS 3 — — which is of order S in the quantity t a + 5S 3 — — . Secondly, as a 

ds a ' ds a 

consequence of these equalities and the following estimates (obtained from (3.5) and (4.1)): 



(5.2) 



dV 

ds a 



Rt c 



(L=(f2))3 



< c 



oS 3 - — n 

ds a 

1 8U s v 



(L2(0))3 



< C 



5 dS 3 
dUsv 



(L 2 (f2))3 



< C 



dSr 



\\di S t(V x v,SO(3))\\ L 2 {Qs) 
\\dist(y x v,SO(3))\\ L 2 {as) 

51/2 

\\di3t{V x v,SO{3))\\ L 2 (as) 
\\dist(V x v,SO(3))\\ L 2 {Qs) 



S 1 / 2 



5n t7 

we deduce that in the quantity Hs((V x v — R) T R + R r (V x v — R)) we can neglect the terms — . 

Now, if in the Green-St Venant's strain tensor of v we carry out the simplifications mentioned above, 
we are brought to replace 

^U g ((V x v) T V x v-I 3 ) by (ti 1 1 2 I n)" r n 4 (£?-(«)) (ti| tain)" 1 
l{(V x v) T V x v-I 3 ) by (t! 1 1 2 I n)-^ 3 ^)^! 1 1 2 | n)" 1 



or 



where the symmetric matrix E s (v) <G (L 2 (£ls)) 3x3 is equal to 



(a 3 Tu(R)+Z u s 3 ri 2 (R) + Zi 2 l~R T p- ■ ti + \z 31 \ 

2 ds 3 2 



(5.3) 



E s {v) = 

r Q/3 (R) 

Zuj3 



S 3 r 2 2(R)+2 2 2 -R T ^— • t 2 + -Z 32 



dv 



\ 



R - — n 



<9R „ dR „ ■ 

-— n • Rt^ + n Rt Q 

ds a dsp 

( dV \ / dV > 

W " RtQ J ' Rt ^ + la^ ~ Rtfi 1 ' Rta 



ds 3 



<7 dV n 



where (ti 1 12 | n) denotes the 3x3 matrix with first column ti, second column t 2 and third column n and 
where (ti 1 12 | n)~ T = ((ti 1 12 | n) _1 ) T . Let us notice that E s (v) belongs to (L 2 (il§)) 3x3 for any deformation 

v(s) = V(si,S2) + S3R(si,S2)n(si,S2) +v(s), for a.c. seftj 

where V G (H 1 ^)) 3 , R G H x (u\ 50(3)) and v G L 2 {u: H^-S, S)) 3 . 
Remark 5.1. From the last estimate in (3.5) we deduce that 



and then we get that the set 



|n 5 (V x w - R)|| (L 2 (n)) 3x3 < CS 



{sefi| |||n 5 (V x w-R)(s)||| > 1} 



□ 



has a measure less than CS . It follows that the measure of the set 

{s g o | dct (n a [Vsu] (s)) < o} 

tends to as 5 goes to 0. 

Now, we introduce the following closed subset B s of (H 1 ^)) 3 x (F 1 (w)) 3x3 x (L 2 (cj; ff^-tf, 5))) 3 
% = {v = (V,R,U) G {H\u)) 3 x (H 1 (w)) 3x3 x (L 2 (lj;H\-5,S))) 3 \ 
R(«i, s 2 ) G 5(9(3), y U(si, s 2 , s 3 )rfs 3 = 0, 

J S3v(si,S2,S3) ■ t a (si,S2)ds3 = 0, for a.e. (si , s 2 ) € w, a =1,2.} 

The last condition on TJ in O5 is not satisfied in general (if v is the warping introduced in Theorem 3.1), 
loosely speaking this new condition will allow to decouple the estimates of v and Z ia (see the proof of 
Proposition 5.2). 

For any veij, we consider v defined by 



(5.4) 



v(s) = V(si, s 2 ) + S3R(si, s 2 )n(si, s 2 ) + v(s), for a.e. s G Cls- 



The deformation u belongs to (L 2 (uj; H 1 (—S,S))) 3 so that, in general, the Green-St Venant's tensor of v is 
not defined. Nevertheless, the tensor field E s (v) belongs to (L 2 (O a )) 3x3 and we set 



(5.5) 



E(v)=E s (v), 



E(v) G (L 2 (n s ) 



3x3 



dV 



Let us point out that if a triplet v satisfies the limit kinematic condition — — = Rt„, then it is easy to 

os a 



obtain 



^ll W ll(L2(n 5 )3 + 



dv 



< \\ E M\\(L 2 (n 5 )) 3 ^ 3 , 



<9R 



ds n 



< 



C 



(L=( W ))3X3 



^I|£(V)||(L2(Q,))3X3 



which permits with some boundary conditions to control the product norm of v in term of | l-^Cv) 1 1 (i 2 (r2 5 )) 3x3 
and S. In order to define an energy which have this property for any v G O5, we are led to add two 
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penalization terms, which vanish as 5 — > 0, to ||-S(v)||^ i2 (n 4 ))3x3- This is why for every deformation v e 
we set 



(5.6) £ 5 (v) = \\E(vW L2 



2 -4- A 3 



dR <9R 

— t 2 — "71 t»l 

OSl OS2 



(i"(u,))3 



- — • Rt 2 - — • Rti 

asi c>s 2 



Proposition 5.2. TTiere exists a positive constant C which does not depend on S such that for all v e 



^IMI(L'(U 4 )S + 

on 2 j_ 

(L 2 (w)) 3 >< 3 (5 2 



<9v 



(L2 (a ,)3 



<9s 



ds 3 

Rt Q 



< I|- B ( V )II(L 2 (0 5 ))3X3 



<9s n 



< S^(v). 

(L 2 (w))3 



Proof. First of all there exists a positive constant C independent of 5 such that 

<9R 



(5.7) 



dR „ 
— n Rt! 

OSi 



i 2 M 



+ 5 



5R „ OR n 
— n-Rt 2 + — n-Rt! 

osi as 2 

C 



£ 2 M 



+ 5 



<9s; 



n Rt 5 



£ 2 M 



+ I|2ii||l 2 (w) + ||2i 2 ||l 2 ( W ) + ||2 2 2||l 2 (w) < ^Y7 2 "ll- E ( v )H(- L2 (^)) 3x3 ' 



rp OV 

R q — • t a + 

ds 3 



i 2 (n a ) 



+ 



<9U 
ds 3 



L 2 (f2 5 ) 



< C||£;(v)|| (L 2 (n4)) 3X8. 



We use the definition of D«5 to estimate the held R T w • t a . Introducing the function R T v ■ t a + s 3 Z 3a , using 
Poincare-Wirtinger's inequality and the first condition on R T v in give 



(5.8) 



R V ■ t a + S^Zsa 



+ 


R T v ■ n 


< 






L 2 (f2 5 ) 



< C5\\E{v)\\ (L 2 (as ))sxs. 



Now we use the second condition on v ■ t a (in the definition of ID>$) in the above estimates and again (5.7) to 

i • t-»T 

get the estimates on R — — and Z 3a 

os 3 



s{| rT k-*-IL™ w + ' i/,|I *- ii ^>} + I rT S-" 



a=l 



(i 2 (^)) 3 



<9v 
<9s 3 



< Cp(v)|| (L2(a4)) 3 X 3. 



Finally (5.8) gives the L 2 estimate on Let us notice that due to the last condition on v in Dj, we obtain 

the same estimates that in the case where v satisfies the limit kinematic condition — — = Rt a . 

os a 

There exist two antisymmetric matrices Ai and A 2 in (L 2 (w)) 3x3 such that 

9R OR 

— = RAi — = RA 2 . 

OSi OS 2 



From (5.7) we get 
II Ain • ti | 



Ain • t 2 + A 2 n • tj J2( ^ + || A 2 n • t 2 \\ L2 



C 



lL 2 (w) 1 ll^ 1 " ' 1 II L 2 (ui) 

Besides there exists a positive constant such 

||Al||(£,2( w ))3X3 + ||A 2 ||(I,2( W ))3X3 



(UJ) - ^72ll S ( V )H(i 2 (^))3x3. 



<c|||Ain-ti|| i2(w) + ||Ain-t 2 + A 2 n-ti| 



r 2 , . + A 2 n • t 2 2 , . 

L (w) H ^IIL 2 (w) 



|Alt 2 — A 2 tl||(I,2( w ))3 j. 



Hence we get 
<9R 



ds i 



(£ 2 («)) 3 



+ 



<9s 2 



(L 2 (w)) 3 >< 3 



< c 



|(L 2 (n 5 ))3X3 + 



<9R t dK t 

dsi 2 ds 2 1 



(L 2 M) 



Due to the estimates concerning the 2j Q and the definition of £a(v) wc finally obtain 



ds a 



We define now the set of the admissible triplets 

jv= (V,R,tJ) eD { | V = 0, R = I 3 on 7o }. 



S7o 



Notice that the triplet 1^ = (</>, I 3 , 0) belongs to H>,5, 70 and it is associated to the deformation v = Id- 

In some sense, the following corollary gives two Korn's type inequalities on the set B,5, 7o with respect 
to the quantity £g(v), the more accurate of which depending on the order of £$(v). 

Corollary 5.3. There exists a positive constant C which does not depend on S such that for all v € O,$, 7o 

IIV- + ||R-I 3 ||^ 1M) 3x3 < J&(v), 
l|V-0||fci (u)) s<c(l + ^ 4 (v)). 

Proof. Recall that R = I 3 and V = <p on 7o, then from Proposition 5.1 we obtain 

C 



I R - I 3|l(ffi( w ))3x3 ^ ^3 



Using the above estimate and again Proposition 5.1 we obtain the first estimate on V — <j> (recall that 



t a = — — ). To obtain the second estimate on V — <j>, notice that ||R — I 3 | [f L 2/ u \\3x3 
os a 

6. Elastic shells 



< C. 



□ 



In this section we consider a shell made of an elastic material. Its thickness 25 is fixed and belongs to 
]0, 2#o]- The local energy W : S 3 — ► M+ is a continuous function of symmetric matrices which satisfies 
the following assumptions which are similar to those adopted in [14], [15] and [16] (the reader is also referred 
to [8] for general introduction to elasticity) 

(6.1) 3c>0 such that \/E G S 3 W {E) > c\\\E\\\ 2 , 

(6.2) Ve>0, 39 > 0, such that V£eS 3 \\\E\\\ < 6 =^ \W(E) - Q(E)\ < e\\\E\\\ 2 , 

where Q is a positive quadratic form defined on the set of 3 x 3 symmetric matrices. Remark that Q satisfies 
(6.1) with the same constant c. 

Still following [8], for any 3x3 matrix F, we set 



(6.3) 



W(F) 



W(^(F T F-I 3 )j if dct(F)>0 
+ oo if det(F) < 0. 
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Remark that due to (6.1), (6.3) and to the inequality |||F T F - I 3 ||| > dist(F, SO (3)) if det(F) > 0, we 
have for any 3x3 matrix F 



(6.4) W(F) > -dist{F,SO(3) f. 

Remark 6.1. As a classical example of a local elastic energy satisfying the above assumptions, we mention 
the following St Venant-Kirchhoff's law (see [8]) for which 

W{F)= i^r(F T F-h)) 2 + ^r((F^F-I 3 f) if det(F) > 
[ + oo if dct(F) < 0. 

In order to take into account the boundary condition on the admissible deformations we introduce the 
space 

(6.5) Vs = {ve(H 1 (Q s )) 3 \v = I d on T . s }. 

Let k > 1. Now we assume that the shell is submitted to applied body forces f K .s € (L 2 (fls)) 3 and we 
define the total energy J Kl s(v)* over JJs by 

(6.6) JM= I W(V x v)(x)dx- f f Kj5 (x) ■ (v(x) - I d (x))dx. 

To introduce the scaling on f K ,s, let us consider / and g in (L 2 (oj)) 3 and assume that the force f Ki s is given 

by 

(6.7) f K , s (x) = S K ' f(s 1 ,s 2 ) + 8 K '- 2 s 3 g(s 1 ,s 2 ) for a.c. x = $(s) G Qs- 
where 



(6.8) k 



2k - 2 if 1 < k < 2, 
k if k > 2. 



Notice that J K ,s{Id) = 0. So, in order to minimize J Kt s we only need to consider deformations v of Us 
such that J Ki s(v) < 0. 

Now from (6.1), (6.3), (6.4), the two Korn's type inequalities (3.1)-(3.2), the assumption (6.7) of the 
body forces and the definition (6.8) of n , we obtain the following bound for ||dist(Va;U, -S'0(3))||l 2 (q 4 ) 

(6.9) | |dist(V x «,SO(3)) | \ L 2 {Qs) <C8 K - X I 2 and / f Ki $ ■ (v — Id) < C5 2k ~ 1 
which in turn imply that 

(6.10) cS 2K - 1 < J K . s (v) < 0. 



* For later convenience, we have added the term / f K ,s(x) ■ Id(x)dx to the usual standard energy, indeed 
this does not affect the minimizing problem for J Kt s- 
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Again from (6.3)-(6.4) and the estimates (6.9) we deduce 

-\\{V x v) T V x v - I 3 ||^ (2s)) 3x3 < J K ,s(v) + f f K ,s ■ (v - I d ) < CS 2K -\ 



'Qs 

Hence, the following estimate of the Green-St Venant's tensor: 



\\{(V x v) T V x v-l 3 }\\ (L2(Q , )3x3 <C6«- 1 / 2 . 



2 \ V *^ / v ^ ||(I,2(Q 4 ))3X3 

We deduce from the above inequality that v € (W 1 ' 4 (Qs)) 3 with 

l|Vx«||(i,4(Q 5 ))3 X 3 < CSi. 

We set 



As a consequence of (6.10) we have 



m K;(5 = inf J K .s(v). 
vev s 



In general, a minimizer of J K ,g does not exist on XJg. In what follows, we replace the elastic functional 
v i — > J K ,s{v) on \Jg by a simplified functional defined on O5 which admits a minimum. 

From now on we assume k > 1. 

7. The simplified elastic model for shells 

The aim of this section is to define a functional J* s on the set B>i, 7o , which will appear as a simplification 
of the total energy J K _g defined on the set \Jg. In order to perform this task, we use the results of Section 5 
and we proceed in three steps. Let us first consider an admissible deformation v satisfying (3.3), decomposed 
as in (3.4) and such that J Ky g(v) < 0. It is convenient to express the energy J K ,s(v) over the domain fig 



(7.1) 



J K A V ) = ( w(l((V x v) T V x v-I 3 ))dct(ti + s 3 ^\t2 + s 3 ^\n)dsids 2 ds 3 

Jq s \Z / OSi OS2 

-J (S K f + S K ~ 2 s 3 g) ■ (v - I d ) det (ti + s 3 J^-|t 2 + s 3 ^-\n)d Sl ds2ds 3 . 



The triplet associated to v by the decomposition (3.4) is denoted v = (V, R, v). The following estimate 
has been proved in Section 6 



< cs K - 1/2 . 

(L 2 (n 6 ))3x3 



Then, for all 8 > 0, the set A e s = {s E ft; \\\Ug((W x v s ) T W x vg - I 3 )(s)||| > 0} has a measure satisfying 



meas(Ag) < C- 



6 2 



Now, according to assumptions (6.2) and k > 1 and the above estimate, in the first term of the total energy 
Jk,s{v) we replace the quantity W (^-((V x v) T V X v — I 3 ^j by q(^-((V x v) T V x v — 1 3 )^ . Following the analysis 

of Section 5, we then replace Q^-((y x v) T V x v — I 3 )^ by Q((ti 1 1 2 | n) _T _E(v)(t 1 1 1 2 In) -1 ) where E(v) is 
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defined by (5.3) and (5.5). At last, we replace dct (ti + S3 — — |t2 + S3 — — |n) by dct (ti |t2 |n) . Setting for 
all 3 x 3 symmetric matrix F 

(7.2) W S (F) = Q ((ti 1 1 2 I n)- T F (ti | t 2 | n)" 1 ) 

all the above considerations lead us to replace the first term in the right hand side of (7.1) by 



(7.3) 



W s (E(v)) dct(ti|t 2 |n)dsids 2 ds3. 



Observe now the term involving the forces in (7.1). We have 



dn 



dn 



(S K f + S K 2 s 3 g) ■ (y - I d ) det (ti + s 3 - — |t 2 + S 3 — \n) d Sl ds 2 ds3 



-28 



K +1 



dsi ds2 
f ■ (V - 4>) + ig • (R - I 3 )n] det(ti|t 2 |n)dsids2 



V-/^.(V-0)[det(^|t 2 |n) + det(t 1 |^|n 



ds\ds2 



<C5 K +2 (||/||(L2 (w))3 + ||fl , ||(I,2( w ))3)(||V - <t>\\(Li{w)Y + ||R- I3||(L 2 ( W ))3X3 + ^72-|kll(L 2 (O s )) 3 )- 

Then, in view of the first estimate in (3.5) we replace the term involving the forces by 

(7.4) £ M (V,R) = <^' +1 £(V,R) 

where 

£(V,R) = 2 f [/.(V-0) + ^-(R-I 3 )nl det(ti|t 2 |n)dsids 2 

+ \j 9- (V-<£)[det (J^|t2 I n) +det (ti|J^ | nj\ds v l^. 



At the end of this first step, we obtain a simplified energy for a deformation weUj which satisfies (3.3) and 
J K ,s(v) < 



J S K 7 P \v)= [ W s (E(-v))det(t 1 \t 2 \n)ds 1 ds 2 ds 3 
Jn s 



3 -(5 K+i £(V,R). 



Indeed the energy J^ pl (v) can be seen as a functional of v defined over D,5 i7o since we have already notice 
that E(y) belongs to (L 2 (Qs)) 3x3 . As a consequence, in a second step we are in a position to extend the 
above energy to the whole set TS>s,-yo an d to P ut 



Vv e 



-CrV) = / W s (E(v)) det(t 1 |t 2 |n)rfs 1 rfs 2 rfs3 - ^' +1 £(V, R). 



As observed in Section 5, the functional J^'™ pl is not coercive on H>8 no . In a third step, in view of Proposition 
5.2 and in order to obtain the coerciveness of the simplified energy, the two terms S 3 

are added to J^™ Pl - 



0R 


dR t 


2 


dsi 2 


ds 2 1 


(£ 2 M) 3 



Rt 2 - — • Rti 

OSi OS2 



L'{ U ) 
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Using all the above considerations, we are able to define the simplified elastic energy on !B>,s, 7o by setting 
for any v in B,5 j7o 



(7.5) 



J S K s(v)= / Vr(£(v))det(t 1 |t 2 |n)ds 1 ds 2 cfo3 + 5 3 



— t -— t " 

dsi ds2 



(L=M)3 



— — • Rt 2 — q — • Rti 

osi os 2 



i 2 M 



S R +1 £(V,R). 



The end of this section is dedicated to show that the functional J* s admits a minimizer on B,5 j70 ■ Let 



v be in B,5 i7o we have 



(7.6) 



£(V,R) < C(||/||(Z,2(u;)) 3 + llffll(L 2 (cu)) 3 )(H V _ 0ll(£ 2 (u))) 3 + ||R- I3||(I,2( W ))3X3). 



The quadratic form Q being positive, the definition (5.6) of £a(v) and (7.5)-(7.6) give 



C£ s (v) - C6 K +1 (||/|| (L 2 (W)) 3 + ||0||( L2(fc))) s)(||V - 4>\\(lw + MR-- l3||(^ (w ))sxs) < J*, s (v) 
Now thanks to Corollary 5.3 and (6.8), we get, if J* )5 (v) < 0( = J* i5 (I<i)) 
(7-7) £,(v) < CT 2K - 1 (||/|| (i2(w)) 3 + || 3 || (LW ) 2 . 



Hence, there exists a constant c which does not depend on 5 such that for any v € B,5, 70 satisfying 5 (v) < 0, 
we have 

cS 2 *- 1 < J' KtS (v). 

We set 



(7.8) 



K,s= jnf ^,«(v). 



As a consequence of the above inequality, we have 



c < 



m 



6 2k- 



< 0. 



In the following theorem we prove that for k and 5 fixed the minimization problem (7.8) has at least a 
solution. 

Theorem 7.1. There exists v$ e B^ 7o such that 



(7.9) 



m 8 K ,6 = jS K Avs) = min J* 5 (v) 



Proof. Since J* = 0, we can consider a minimizing sequence v„ in B,5 j7o such that J^s( v n) < and 



From (7.7) we get 



n— J-+00 ' 



^«(V„) < C(5 2K (||/(L2( w ))3 + ||5|| (L 2 (s:2)) 3) 4 . 
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Thanks to Corollary 5.3 and Proposition 5.2, the above estimate show that there exists a subsequence still 
denoted n such that (recall that ||Rn||(L°°(a;)) 3 x 3 = \/3) 

V n ^V s weakly in (H 1 ^)) 3 
R„ Kg weakly in (fl" 1 (w)) 3x3 
R„ — > Rg strongly in (L 2 (uj)) 3x3 and a.e. in u 

v n ^v s weakly in (L 2 (u>; H 1 (—6, 6))) 3 . 

Then setting v 5 = (Vg, R.g, vg) e %, 70 , we get 

E(v n ) -± E(vg) weakly in (L 2 (f! 5 )) 3x3 , 

^ ■ R„t 2 - ^ ■ R„t! - ^1 . R , t2 - p. . R, tl weakly in (L 2 {u)) 3 . 

OSi OS 2 OSi OS 2 

Now, passing to the limit inf in J* g (v n ), we obtain 

m s Kt g < J^sivg) < lim J*, 4 (v„) = lim J^(v n ) = ro' tf . 



8. Asymptotic behavior of the simplified model. Case k = 2. 

In this section we study the asymptotic behavior of the sequence (yg) of minimizcr given in Theorem 7.1 

and we characterize the limit of the minima — ^— as a minimum of a new functional. AS usual, to perform 

<r 

this task, we work on the fixed domain il and we use the operator 11,5 defined in Section 4. We denote D the 
following closed subset of ID>i l70 (i.e. Hs,-y for <5 = 1 or B>i j70 = 11,5 (B^)): 



{v = (V,R,V)e©i, 70 I |£=Rt Q }. 



Notice that V € (H 2 (uj)) 3 . Then we define the following functional over D 

(8.1) J 2 (v) = J Q((ti 1 1 2 | n)- T E(v)(t! 1 1 2 | n)- 1 ) det(ti|t 2 |n) - C(V, R). 

where 



(8.2) B(v) 



Ss |£n.R t2 iR^.t^ 

„ <9R M 1„ T <9F 

* S 3 — n-Rt 2 -R T — -t 2 

<9s 2 2 9S3 



v * * R ^- n y 



dV 



As in Theorem 7.1 we easily prove that there exists v 2 = (V 2 , R 2 , V2) € D such that 
(8.3) m s 2 = J 2 (v 2 ) = minj 2 (v). 

Theorem 8.1. We have 

m 2 s 

to 2 = lim — -j-. 
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Moreover, let v$ = (Vs,Rs,vs) € D^ 7o be a minimizer of the functional J 26 {-), there exists a subsequence 
still denoted S such that 



(8.4) 



I <5 2 



V s — > V strongly in (H 1 (w)) 3 , 
R«5 — > R strongly in (H 1 (lu)) 3x3 , 

-zZififi — > strongly in L 2 (u>), 

Rs(vs) — > V Q strongly in (L 2 (uj; if^-l, l))) 3 . 



The triplet v = (Vo,Ro, Vo) belongs to D and we have 

m s 2 = ^;(vo). 

Proof. For all v = (V, R, V) e D, we have (V, R, v s ) € D,5, 7o where 

vs(si, s 2 , S3) = S 2 V(si, s 2 , for a.c. (si, s 2 , S3) G f2j. 

<9R <9R 

Using the fact that v e B, which implies that — — n • Rt 2 = — n • Rti, we have 

OSi Os 2 

(8.5) J 2,a(V,R,"i) = J Q^ tl 1 t 2 I n)- T E(v)(t! 1 1 2 | n)" 1 ) det(ti|t 2 |n) - £(V, R) = J 2 (v). 

m 2 <5 

Then, taking the minimum in the right hand side w.r.t. v € B, we immediately deduce that ' < m 2 . 
We recall that vs = (Vs,Rs,vs) € Ds ll0 is a minimizer of J| 

z ™2,<5 ^(v*) J 2 S ,5( V ) 

c - ~TT = Fi = mm — 73 — 



and moreover with (7.7) 



£<$(v{) < C , (5 3 (||/( i 2( w)) 3 + ||5||(x,2(n))3) . 



Thanks to the estimates in Proposition 5.2, Corollary 5.3 and the above estimate we can extract a subsequence 
still denoted S such that 

Vs^Vo strongly in (H 1 ^)) 3 , 
Rs — 1 Rq weakly in (iJ 1 (w)) 3x3 and a.e. in u>, 



(8.6) 



^n 5 (iJ 5 ) - Vo weakly in (L 2 ( W ; tf^-l, l))) 3 , 



1 



Zia,o weakly in L 2 (u>) 
■ OR, . 0R S 



\0 Sl 2 



l/OVs 



{ 5 V 9si 



R 5 t 2 - 



0s 2 
0s 2 



ti) y weakly in (L 2 (cj)) 3 , 
• R<sti^ — >■ X weakly in L 2 {lo). 



0V a 



Then from the fifth convergence we obtain = Rot a . So we have Vo € (H 2 (lj)) 3 and v = (V ,Ro, Vo) 
belongs to D. From the above convergences, and upon extracting another subsequence, we also get 

in 5 (£( V(5 )) - E weakly in (L 2 ^)) 3 * 3 
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where 



E — 



£3 — " n ■ Roti + Zufi S3 " n ■ Rot2 + £12,0 ttR-o 



9si 



9Ro 
<9si 

<9s 2 



V 



9Sj 

S , 3^7^n • R t 2 + Z 2 2,o 2"^° "95^ 

R ° asT' 



•ti 



n 



with 



Wo — Vq + Sa^i^Roti + 5*32,32, oRot 2 
Due to the expression of J| 5 we have 



<5 3 



= J g((ti 1 1 2 I n)- T ln,(£(v«)) (ti 1 1 2 I n)- 1 ) det(ti|t 2 |n)d Sl ds 2 rf53 

£(V 5 ,R 5 ). 



1 

+ S2 



&Rs. 


ORs t 


2 


OSi 


ds 2 


(£ 2 M)3 



— — • R«5t 2 — — • R«5ti 

OSi os 2 



L=»(u») 



With the convergences (8.6), since Q is quadratic and thanks to the expression of £, we are in a position to 
pass to the limit-inf in the above equality which gives 

J Q((ti 1 1 2 I n)- T E (t! 1 1 2 I n)- 1 ) det(ti|t 2 |n)d*ids 2 dS3 + 11*11*^) + IMlfiW - £(Vo, Ro) 

< urn — — = urn 



(5^0 



J 3 



(5^0 



<5 3 ' 



Hence we get 



/ Q((t 1 |t 2 |n)- T Eo(t 1 |t2|n)- 1 )det(ti|t 2 |n)dsids 2 rf5 , 3 -£(V ,Ro) < Mm 



2,5 



<5 3 



First, notice that if v = (V, R, V) € B then V satisfies 
"1 »T7 



tttt - ( s i ) s 2 1 S3) ■ t a (si, s 2 )(Sf - l)dS 3 = for a.e. (si, s 2 ) G w. 
1 0^3 

Now we apply Lemma A with a = ( ^— -n ■ R ti, ^— • Rot 2 , ^r-^n 1 Rot2~) , b = (Z\\ , Z\ 2 0, -Z 22 ), 

V asi asi os 2 / 



c = (7-R0 ' ^R-o ' t 2 ,Ro -773— ' n ) and with the quadratic form defined by 

V 2 003 2 003 063 / 



Q m (a,b,c)= J Q((ti|t2|n)- T Eo(ti|t2|n)- 1 )det(ti|t2|n)dS 3 for a.e. («i,s 2 ) 



We obtain 



(8.7) minj 2 (v)< / Q((ti 1 1 2 | n)- T E (ti 1 1 2 | n)" 1 ) dct(t 1 \t 2 \n)d Sl ds 2 dS 3 - £(V , Ro). 



in: 



Hence m 2 < lim ^ . Recall that we have —4— < m%, so we get 



m 



2,5 



TO 2 5 

lim — -f- = m s 2 . 
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Finally, from convergences (8.6) we obtain Z ia , = 0, X = Y = and moreover we have the strong 
convergences in (8.4). □ 

9. Justification of the simplified model. Case k = 2 

In this section, the introduction of the simplified energy is justified in the sense that we prove that both 
the minima of the clastic energy and of the simplified energy have the same limit as 5 tends to 0. 
Theorem 9.1. We have 



lim 



m 2 ,6 



lim 



'2,5 



to;. 



(5^0 S 3 <5^0 S 3 

Proof. 

Step 1. In this step we prove that m% < lim — . Let (vs)„, x ^ x be a minimizing sequence of deformations 

S^q 3 'U<<>S<>0 

belonging to XJs and such that 



(9.1) 



From the estimates of Section 6 we get 



lixn = lim M v <) 



5-s-O 



S 3 s^o S 3 



(9.2) 



f ||dist(V^ 4 ,SO(3))|| ia(e4) <C6 3 ' 2 , 

\ \\l{VxvJV x v s - I S }\\ (LHQ5))3X3 < OS 3 ' 2 , 

{ ||Va;Ui||(I,4(Q 4 ))3X3 < C^^. 



We 



1 f 5 

still denote by Vs(si,S2) = — ~ / v$(si, s 2 , s^ds^ the mean of v$ over the fibers of the shell. Upon 

2d Js 

extracting a subsequence (still indexed by 6), the results of [6] show that there exist V G (H 2 (lj)) 3 , R G 
(tfV)) 3 * 3 withR(si,s 2 ) G 50(3) for a.e. (si,s 2 ) G w, Z Q) g e L 2 (w) andFe (L 2 (w; ff^-l, l)) 3 satisfying 



V(si, s 2 , S 3 )dS 3 = for a.e. (si,s 2 )6w, ^ — = Rt c 



,9.3, /_; 

together with the boundaries conditions V = 0, R = I 3 on 70, and with the following convergences 
'Rj(» 4 )— >-V strongly in (fT^fi)) 3 , 

n 5 (v x w5) — 
n 5 ( U5 - v s ) 



(9.4) 



where 



ILj(V x w 5 ) — > R strongly in (L 2 (Q)) 3x3 , 
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S 3 (R-I 3 )n strongly in (i 2 (0)) 3 , 
ILj((V x w 5 ) T V x t; 5 - I 3 ) ->■ (ti 1 1 2 I n)- T E (ti 1 1 2 I n)- 1 weakly in (L 2 (^)) 9 , 



' S 3 — n-Rti+Zn S 3 — n-Rt 2 + Z 12 -R — • ti 

asi asi 2 aA 3 



E 



n <9R „ „ 1„ T 5V 
S 3 — n-Rt 2 + Z 22 -R T — -t 2 

as 2 2 9S3 



Now, recall that 

■h,s{ v s) 



(9.5) 



6 3 



^^w(in 5 ((v x ^) T v x ^-i 3 ))n 5 (det(v$))-^ ^ f K , s • (v s - i d ) 
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In order to pass to the lim-inf in (9.5) we first notice that dct(V$) = det(ti|t 2 |n) + S3 dct ^J^-|t 2 | nj + 

s^det (ti|— — | n) + S3 dct ( — — | — — | n) so that indeed 11$ ( det(V$)) strongly converges to det(ti|t 2 |n) 

V 0S2 ' \os\ 0S2 ' 

in L°°(0) as 8 tends to 0. 

We now consider the first term of the right hand side. Let e > be fixed. Due to (6.2), there exists 
9 > such that 

(9.6) V£ e S 3 , |||£7|||<0, W(E) > Q(E) - e\\\E\\\ 2 . 

We now use a similar argument given in [5]. Let us denote by xi the characteristic function of the set 
A% = {s e 0; \\\U s ((\7 x v s ) T \7 x v s - h)(s)\\\ > 0}. Due to (9.2), we have 

(9.7) meas(Ag) < C^. 
Using the positive character of W, (9.2) and (9.6) give 

J ^w(n 5 (v x v 5 ))\n 5 (dct(v<i>))\ > ^w(^n s ((v x v 5 ) T v x v 5 - i 3 ))(i - x*)n 5 (det(v$)) 

> J^Q(^n s ((W x v s ) T W x vs - I 3 )(l - x?))n 4 (det(V$)) - Ce 

In view of (9.7), the function \s converges a.e. to as S tends to while the weak limit of -^\s{(y xVs) T ^ X V8— 
I3XI — Xs) 1S given by (9.4). As a consequence and also using the convergence of II,5(det(V<I>)) obtained 
above, we have 

lim / l|?(n 5 (V x ^))n 5 (det(V$)) > / Q((t 1 |t 2 |n)- T E(t 1 |t 2 |n)- 1 )det(ti|t 2 |n)-C£. 
s^oJn Jn v ' 

As e is arbitrary, this gives 

(9.8) lim f ^rW(n 5 (V x v 5 ))ll s (det(V$)) > [ Q((ti | t 2 | n)~ T E (ti 1 1 2 1 n) -1 ^ det(ti|t 2 |n). 
Using the convergences (9.4), it follows that 

where £(•, •) is defined by (8.5). From (9.5), (9.8) and the above limit, we conclude that 

(9.9) lim ^ = lim > f Q ( {tl | t 2 | n)~ T E (ti 1 1 2 | n)" 1 ) det(ti|t 2 |n) - £(V, R). 
Proceeding as in the proof of (8.7) in Section 8, we get 

/ Qf(ti 1 1 2 I n)- T E (ti 1 1 2 I n)- 1 ) det(ti|t 2 |n) - £(V, R) > mmj 2 (v) 



TO 9 . 



Finally we have proved that m s 2 < lim m ^' 5 . 

8^0 <5 d 
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Step 2. In this step we prove that m 2 > lim 

Let us now consider a minimizer vo = (VojR-Oi^o) € ID> of Ji and the sequence ((V 4 ,H5,V 4 )) Of 
approximation of v given by Lemma C constructed in the Appendix. The deformation v$ is now defined by 

(9.10) vg(s) = Vs(si,s 2 ) + s 3 R s (s 1 ,s 2 )n(s ll s 2 ) + 6 2 Vg(si,S2, y), for s e £ls- 

Step 2.1. Estimate on \\ILs(V x vs - R5)||(L~(n))3x3 and | \dist(\7 x v s , 50(3)) 
From (9.10) and trough simple calculations, we first have 



(9.11) 



then 



,_ _ . dVs „ 9R,5 c2 dVs ,_ _ . 9n 

(V^ - R$)t Q = R<5t Q + s 3 — — n + 6 — (\/ x v s - R<5) s 3 ^ — 

os a as a as a os a 

(Vxvs - R«5)n = 5^^-, 
003 



U 5 (W x v 5 - R 5 ) ■ n 4 (v s $) 

= R^ti+S^- — n + — R,5t 2 + S^- — n + <5 2 - — <S— - 

V asi osq, os q c*s 2 c*s Q c*s Q 0^3 

Thanks to (2.3) and the estimates of Lemma C in Appendix we obtain 

(9.13) ||ILj(V x ^ - R*)||(ioo ( n))3*s < \ 
and we deduce that there exists a positive constant Cq such that 

(9.14) \\H s ((V x v s ) T V x v s -I 3 )||(Loc ( n))3x3 < C . 
Again using the estimates in Lemma C we get 

\\dist(V x v s ,SO(3))\\ L oo (u) < l - 

and then we obtain 

(9.15) for a.c. s e ft 5 det (V x v 5 (s)) > 0. 

Step 2.2. Strong limit of ^IL; ((V x v s ) T V x v s - I 3 ). 

Thanks to the estimates and convergences of Lemma C and (9.12) we have 

(9.16) ||n 5 (V x w 4 - R«)||(i,2(n))3x3 < CS. 

We write the identity {V x v s ) T V x v s - I 3 = {V x v s - R s ) T Rs + Rj (V x v s - R s ) + {V x v s - R s ) T {V x v s 
R s ) + (R s - R) T R<5 + R T (R,5 - R). So, from (9.13) and (9.16) we get 



(9.17) 



\\n-s((Vxvs) T v x v 5 - 13) I l(i 2 (n)) 3 >< 3 < cs. 
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In view of (9.11), the strong convergences of Lemma C and (9.16) we deduce that 



(9.18) 



1 9R 

-ILj((V x ^-R)t a ) -^S 3 —n strongly in (L 2 (0)) 3 

ln s ((V x vs - R)n) — > |K . n strongly in (L 2 (0)) 3 

v. C03 



Now thanks (9.13) and the strong convergences (9.18) we obtain 

4=n«(V x v 5 - R) — > strongly in (L 4 (0)) 3 
Vo 

and then using again Lemma C, (9.18) and the above decomposition of (V x vs) T V x vs — I3, we get 

(9.19) ^((V^fV^-Is) — ^ (t 1 |t 2 |n)-' r E( V o)(ti!t 2 |n)- 1 strongly in (L 2 (ft)) 3x3 , 
where E(vo) is given by (8.2). 

Step 2.3. Let e be a fixed positive constant and let 9 given by (7.2). We denote Xs the characteristic 
function of the set A 9 5 = {s e O; \\\Tl s ((V x v s ) T V x v s - I 3 )(s)||| > 9}. Due to (9.17), we have 

S 2 

(9.20) meas(^) < C-^ 

and from (9.15) we have det (V x vs(s)) > for a. e. s e Clg. Due to (6.2), (6.4) and (9.19) we deduce that 

lim" / l(l- X f)W(n s (V x i; s ))n 5 (det(V$)) < f g((t x 1 1 2 I n)- T E(v )(t 1 1 1 2 | n)" 1 ) det(t 1 |t 2 |n) 

+ e / |||(t 1 |t 2 |n)- T E(v )(t 1 |t 2 |n)- 1 ||| 2 dct(t 1 |t 2 |n) 

where E(v ) is given by (8.2). Notice that there exists a positive constant C\ such that for all E € S3 
satisfying 9 < \\\E\\\ < Co we have 

W(E) < Ci|||£7|||. 

Thanks to (6.3), (6.4), (9.17), the strong convergence (9.19) and the weak convergence ^xt ^ in L 2 {VL) 
we obtain 

lim / 4x?W(n«(Vx««))lItf(det(V*)) < Ci lim / - X l \\\^U s ((V x v s ) T V x v 5 - I 3 )||| II*(det(V$)) = 
Hence for any e > we get 

En I lT?(n 5 (V^))n 5 (det(V$)) < [ Q((t 1 |t 2 |n)- T E(v )(t 1 |t 2 |n)- 1 )det(t 1 |t 2 |n) 

+ e / |||(t 1 |t 2 |n)- T E(v )(t 1 |t 2 |n)- 1 ||| 2 det(t 1 |t 2 |n) 
Jn 

Finally 

(9.21) lim"/ ^-W(n s (V x v s ))n s (det(V<S>)) < [ q((U | t 2 | n)- T E(v )(ti 1 1 2 | n)" 1 ) dct(ti|t 2 |n). 
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As far as the contribution of the applied forces is concerned, we use the convergences of Lemma C to obtain 



(9-22) Jim (1 | q h.s ■ (vs I d j) = C(V, R). 

From (9.21) and (9.22), we conclude that 

< ^Q((ti 1 1 2 | n)- T E(vo)(t! 1 1 2 | n)- 1 ) det(ti|t 2 |n) - £(V,R) = J 2 (v ) = m s 2 . 



ti J2,g( v s) 

hm -= — 

<5^0 5 3 



Then we get lim — ^ < mi,. □ 

10. Alternative formulations of the minima m* >5 and m|. 

In the following theorem we characterize the minimum of the functional J® 5 (-) over IUj 70 , respectively 
^2 over D, as the minima of two functionals which depend on the mid-surface deformation V and on the 
matrix R which gives the rotation of the fibers. 

The first theorem of this section shows that the variable v can be eliminated in the minimization problem 
(7.9). 

We set 

E= {(V,R) e (#V)) 3 x (#V)) 3x3 I V = & R = I 3 on- 
R(si, s 2 ) e SO(3) for a.e. (si,s 2 )ew|. 

We recall (see (5.3)) that for all (V, R) € E we have set 



■7o, 



1 



Rt„ ) Rt +(^-Rt^) -Rt, 



Z?, a — q — • Rn 



r ^ R ^2{^ n - Rt ^ + ^ n - Rt 4- 



Theorem 10.1. Let = (Va, R5, V$) e D5 7o smc/i fftaf to* 5 = J S K 5 (vs) = min J* 5 (v). We have 

veB Srm 



(10.1) 



m%,6 = ^m(V«,R«) = min^ ^^(V.R) 



where 



(10.2) { 



'j-*, s (V,R)=<J 3 / a a/3a , /3 ,r^(R)r a , /3 ,(R)+6 / 



+ <5 3 



9R t 9R t - 

ds 1 9s 2 



(L»(a,))3 



— — • Rt 2 — — • Rti 

OSi OS 2 



i 2 M 



S K +1 £(V,R). 



TTie a Q/3c / p> and b iai > a > are constants which depend only of the quadratic form Q and the vectors (ti, t 2 , n). 
Proof. We have 



rn 



K,5 



mm 

v€Di T , 



6^ Q((t 1 |t 2 |n)- T n 5 (^(v))(t 1 |t 2 |n)- 1 )det(ti|t 2 |n)rf Sl rfs 2 dS3 

-<5 K ' +1 £(V,R) 



<9R 


da 


2 




dV 


— — 1 2 

OSi 


ds2 1 


(£ 2 M) 3 




dsi 



dV 2 
■ Rt 2 — t; — ' Rti 

ds 2 L 2 (u>) 
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In order to eliminate v, we first fix (V, R) G E. We set 

J Q((ti 1 1 2 | n)- T n 5 (^(v))(t! 1 1 2 | n)- 1 ) det(ti|t 2 |n)dS 3 = 3™(a, b, c + d) 



where 



/r n (R)\ /Zu 
a = 6\ r 12 (R) , b = Zi 2 

\r 22 (R)/ \z 22 



z 31 




2 R ^sr'* 1 ^ 

2 R ^T' t2 



T an^(i;) 
R "asT 



/ 



and we apply Lemma B in Appendix to obtain the theorem. 



□ 



The next theorem is similar to Theorem 10.1 for the limit energy and the minimization problem (8.3). 
We set 

E Km = {(V,R)eE I |H-=Rt a , a =1,2}. 
Theorem 10.2. Let v = (Vo, Ro, Vo) G D such that m% = Jii^o) = minj^fv). We have 



(10.3) 

where 
(10.4) 

The a 



12 =^ 2 (V ,Ro) = min J- 2 (V,R) 
(V,R) G E Hm 



^ « aPa 'p' ' R ^) ' R Vj - R) 



^ pi are the same constants as the one in Theorem 10.1. 

Proof. We proceed as in Theorem 10.1. In order to eliminate V, we fix (V, R) G Ejj m and we minimize 
the functional ^(V, R, •) over the space V. Thanks to Lemma B in Appendix we obtain the minimum with 
respect to V and then the new characterization of the minimum to 2 . □ 

Of course, for all (V, R) G E; im , we get 

J-|,i(V,R) =(5 3 J- 2 (V,R). 

Let us give the explicit expression of the limit energies T S K s and J" 2 in the case where S is a developable 
surface such that the parametrization is locally isometric 

V(si,s 2 )guJ ||t Q (si,s 2 )|| 2 = 1 ti(si,s 2 ) • t 2 (si,s 2 ) = 0. 



We consider a St Vcnant-Kirchhoff's law for which we have 



W(F) 
so that Q = W = W s . 



±(tr(F T F~I 3 )) 2 + ^tr((F T F-I 3 ) 2 ) if det(F) > 
+ 00 if dct(F) < 0, 
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Expression of J 7 * s . For any v = (V,R, v) G D«5 l7o , the expression (7.5) gives 



(10.5) 



J s K ,sM = * J n [^N^(v))) 2 + »tr((E(v)f) 
+ 5 



+ 5' 



— — • Rt 2 — q — Rti 

asi ds 2 



<9R dR 2 
o — *2 — t; — ti 

asi as 2 (£ 2 M) 3 

(5 K ' +1 £(V,R). 



where E(v) is defined by (5.3). It follows that the elimination of V in Theorem 10.1 gives the partial 
derivatives of V with respect to S3 



(10.6) 



dv 
ds 3 
dv 
ds 3 

v ^(.,., S3 )-Rn 



(., .,s 3 ) • Rt 2 



£31 



(s 2 + l(4-s 2 )) 
\-j^( s 3 [ r n( R ) + r 22 (R)] + [z u + z 22 



\ 



and then 



t' KiS (v,b)= ^ [ (r Q ,(R)) 2 + Krn(R) + r 22 (R)) 2 " 



a,0=l 
2 



a,/3=l 

dR dR I! 2 

IT" *2--5— *1 + 5 

asi as 2 II(l 2 (w))3 
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dV „ dV 2 



- — • Rt 2 - — • Rti 

dsi ds 2 



L2( W ) 



-<5 K ' +1 £(V,R). 



Expression of J" 2 . For any v = (V, R, V) G D, the expression (8.1) gives 

J2(V) = X [^(ME(v))) 2 +Mtr(( E ( v )) 2 )] -AV,R) 

where E(v) is defined by (8.2). It follows that the elimination of V in Theorem 10.2 is identical to that of 
standard linear elasticity (see [18]) hence we have 



(10.7) 
and then 



V(.,.,S 3 ) = - 



2(1 -f) 



5.? 



rn(R) + r 22 (R) 



Rn 



2 

j- 2 (v,r) = — - / [(1 — v) ]T (r^(R)) 2 + ^(r 11 (R) + r 22 (R)) 2 l -£(v,R). 
^ A v > J " a ,p=i 

Remark 10.1. In the case of a St-Venant-Kirchhoff material a classical energy argument show that if 
( v s) 0<s<So is a sequence such that 

s r -h.s{ v s) 
m 2 = hm — — — , 

then there exists a subsequence and (Vo,Ro) G E, which is a solution of Problem (10.3), such that the 
sequence of the Green-St Venant's deformation tensors satisfies 

^Il s ((V x v s ) T V x v s - I 3 ) — ► (t! 1 1 2 I n)- T E(v )(t 1 1 1 2 I n)- T strongly in (L 2 (f!)) 3x3 , 
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where E(v ) is defined in (8.2) with Vq given by (10.7) (replacing R by R ). 

dV dV 

Remark 10.2. It is well known that the constraint — — = Rti and — — = Rt2 together the boundary 

osi os 2 

conditions are strong limitations on the possible deformation for the limit 2d shell. Actually for a plate or as 
soon as S is a developable surface, the configuration after deformation must also be a developable surface. 
In the general case, it is an open problem to know if the set E;i m contains other deformations than identity 
mapping or very special isometries (as for example symetries). 

Appendix. 



Lemma A. Let Q m be the positive definite quadratic form defined on the space 



x K" x 



(L 2 (-l,l)) 3 by 



V(a,b,c)£l 3 xI 3 x(L 2 (-l,l)) 3 , Q m (a,b,c) = y A(S; 



3) 



^S 3 ai + bi^ 
S3&2 + b2 
S 3 a 3 + b 3 

ci(S 3 ) 
c 2 (S 3 ) 
V c 3 (S 3 ) J 



f 5 3 ai + bi ^ 
S" 3 a 2 + b 2 
S 3 a 3 + b 3 

ci(S 3 ) 
c 2 (S 3 ) 
V c 3 (S 3 ) J 



dS 3 



where A (S3) is a symmetric positive definite 6x6 matrix satisfying 



(Al) 



A(S 3 ) = A(-5 3 ) for a.e. S 3 &}-!,![ 



and moreover there exists a positive constant c such that 



(A.2) 



V£ G K 6 , A(S 3 )£ ■ £ > c|C| 2 for a.e. S 3 e] - 1, 1[. 



For all ael 3 , we have 



min Q m (a, b, c) = min Q m (a, 0, c) 

(b,c)en 3 x(L 2 (-l,l)) 3 cGL 2 



where 



L 2 = {ce(L 2 (-l,l)) 3 | J^c a (S 3 )(S 2 3 -l)dS 3 = 0, ae{l,2}}. 



Proof. We write 



A(S 3 ) = 



Ai(5 3 ) : A 2 (S 3 ) 



.A^(S 3 ) ! A 3 (S 3 ), 



where for a.e. S 3 e] — 1, 1[, Ai(S 3 ) and A 3 (S 3 ) are symmetric positive dchnitc 3x3 matrices. The both 
minimum are obtained with 

c (S 3 ) = -A^(S 3 )A^(S 3 )S 3a , b = 0. 



We have 
(A3) 



Q m (a,0,c ) = ( y Sf (A!(5 3 ) - A 2 (,S 3 )A 3 - 1 ( > S 3 )A^(5 3 ))d5 3 )a • 



a. 



□ 



In the following lemma we use the same notation as in Lemma A. 
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Lemma B. Let a, b be two fixed vectors in R 3 and let d be a fixed vector in R 2 x {0}. We have 



{B.l) min Q m (a,b,c + d) 

c£L 2 



( | 1 i 53 2 [A 1 (^)-A 2 (53)A 3 - 1 (53)A^(5 3 )])a.a + Q m (b,d) 



where Q m is a positive definite quadratic form which depends only on the matrix A. 

Proof. Through solving a simple variational problem, we find that the minimum of the functional c 
Q m (&, b, c + d) over the space L 2 is obtained with 



where e e R 2 x {0} 



c(5 3 ) = d A^(S 3 )A^(S 3 )(S 3 a + b) + (S 2 - l)A^(S 3 )e 



1\ /(T 
e = ( i.e, , 2 e_>. e, = | I , e 2 = I 1 



is the solution of the system 



e a = 0, a = 1, 2. 



'|d- ( |\s 3 2 - l)A,^ 1 (,S3)A 2 n (53)^ 3 )b+ ( |V 3 2 - 1) 2 A 3 - 1 (5 3 )^ 3 ) 

Notice that the matrix J (S 3 — l) 2 A^ 1 (S 3 )dS 3 is a 3 x 3 symmetric positive definite matrix. Replacing c 
and e by their values we obtain (B.l). □ 
Lemma C. Let(V,R,V) beinT> 2 , there exists a sequence ((V 5 ,Ri,Fi)) o/(iy 2 ' 00 (o;)) 3 x(W 1 ' 00 (w)) 3x3 x 



(W 1 ' 00 ^)) 3 sucft toot 



(C.l) 



with 



(C.2) 



Vs = 4>, H ) = T s on 70, = 0, on 7 x]-l,l[, 



Vs — >V strongly in (H 2 (lj)) 3 
Rs — > R strongly in (H 1 (uj)) 3x3 
(Rj-R) — >0 strongly in (L 2 (uj)) 3x3 

RstSj — > strongly in (L 2 (lu)) 3 
V 5 — >V strongly in (L 2 (u; H\(-l, l))) 3 , 
■ strongly in (L 2 (Q)) 3 , 



1 

6 

IfdVs 

8 \ds n 



.dV 5 



ds 



and moreover 



1 

< -. 



(C.3) 



<9s Q 



^(R 5 ,SO(3))|Uoc (w) 
The constant c 1 is given by (2.3). 



1 

< -. 



(L°°(u;))3 
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Proof. For h > small enough, consider a C^(R 2 )-function ip h such that < tph < 1 

J iph(si,s 2 ) = 1 if dist((s!, s 2 ), 7o ) < ft 
\ i>h(si,s 2 ) = if dist((si, s 2 ), 70 ) > 2ft. 

Indeed we can assume that 



C C 
(C.4) ll^/illw 1 . 00 ^ 2 ) < llVVtllw 2 ^^ 2 ) < ^2- 

Since w is bounded with a Lipschitz boundary, we first extend the fields V and R n = Rn into two fields of 
(i? 2 (R 2 )) 3 and (iJ 1 (M 2 )) 3 (and we use the same notations for these extentions). We define the 3x3 matrix 
field R' G (H 1 (M 2 )) 3x3 by the formula 

(C.5) R=(-|-|R n )(t lN „) . 

By construction we have — — = R t Q in R 2 and R = R in to. At least, we introduce below the approxima- 
te* 

tions Vh and R/j of V and R as restrictions to uj of the following fields defined into R 2 : 
V' h {si,s 2 ) = ——J / V(s 1 +t 1 ,s 2 +t 2 )dt 1 dt 2 , 

y^h JB(0,3h) 

R/ l (si,s 2 ) = — i-r / R (si +ii,s 2 +t 2 )dt 1 dt 2 , 

wh £ JB(0,3h) 

and 

(C.7) Vfc = ^fc+Vfc(l-^ fc ), R h = I 3 V^ + Rl(l-^), in w. 

Notice that we have 

V,; g (VK 2 '°°(R 2 )) 3 , R^ g (W 1 ' 00 (R 2 )) 3x3 ) 



(T'.(i) < , r ( °' 3h) a.e. (si, S2 ) £M 2 

i,3/i) 



(C.8) 

V ft G (VK 2 '°») 3 , R ft G (W 1 ''°°(oj)) 3x3 , V h = cf>, R h = I 3 on 7o . 



Due to the definition (C.5) of R and in view of (C.6) we have 
(C.9) 

and thus using estimates (C.4) 



V h — ► V strongly in (ff 2 (R 2 )) 3 , 
R h — > R' strongly in (iJ 1 (M 2 )) 3x3 



f V^V strongly in (ff 2 M) 3 
(CIO) < 

I R ft — > R strongly in (if^u;)) 



3x3 



Moreover using again (C.6) and the fact that R — R^ strongly converges to in (H 1 (M. 2 )) 3x3 we deduce 
that 

^(R.' h - R') — >• strongly in (L 2 (R 2 )) 3x3 
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and then together with (C.4), (C.5), (C.7) and (CIO) we get 



(R/i - R) — > strongly in (L 2 (oj)) 3 *\ 



- R fc t a ) — ► strongly in (L 2 (oj)) 3 



We now turn to the estimate of the distance between R/j(si,s 2 ) and SO (3) for a.e. (si,s 2 ) € w. We apply 
the Poincare-Wirtinger's inequality to the function (1*1,1*2) 1 — ► R (1*1,1*2) in the ball B((si, s 2 ), 3ft). We 
obtain 

/ 1 1 |R (1*1,1*2) -Rl(si,s 2 )||| 2 diiidM 2 < C/\ 2 ||VR'||^ 2(jB((siiS2)i 3 M)) 3 

JB((si,s 2 ),3(i) 

where C is the Poincare-Wirtinger's constant for a ball. Since the open set oj is boundy with a Lipschitz 
boundary, there exists a positive constant c(uj), which depends only on oj, such that 

|(S((si,S2),3ft)\S((si,s 2 ),2ft)) nw| > c(w)ft 2 . 

Setting m,h(si,s 2 ) the essential inhmum of the function (1*1,1*2) l_ |||R(i*i, 1*2) — R- ft (si , S2) 1 1 1 into the set 
(i?((si, 82), 3ft) \ S((si, S2), 2ft)) n oj, we then obtain 

c(w)ft 2 m ft (s 1 ,s 2 ) 2 < Cft 2 ||VR'|| 2 i2(i3((si!;S2)! 3, l)))3 

Hence, thanks to the strong convergence of R h given by (C.9), the above inequality shows that there exists 
h which does not depend on (si, s 2 ) € oj such that for any ft < h 

dist(TL h (si, s 2 ), SO(3)) < 1/8 for any (si, s 2 ) e ZZJ. 

Now, 

• in the case dist((s\, s 2 ), 70) > 2ft, (si, s 2 ) € oj, by definition of R/j and thanks to the above inequality 
we have dist(R h ( Sl , s 2 ), SO(3)) < 1/8, 

• in the case disi((si,s 2 ),7o) < h, (s\,s 2 ) € oj, by definition of R^ we have R/j(si,s 2 ) = I3 and then 
dist(R h {s 1 ,s 2 ),SO(3)) = 0, 

• in the case ft < dist((si, s 2 ), 70) < 2ft, (si,s 2 ) € w, due to the fact that R = I3 onto 70, firstly we 

have 

l|R - I 3|l(L2( W6h . yo ))3X3 < Cft 2 ||VR \\f L 2 {uBhno))3 x3 

where Wfch j7o — {(si,s 2 ) e M 2 I dist((si,s 2 ),7o) < fcft}, fc e N*. Hence 

ll R h - I 3|l(L2( W3h , 70 ))3x3 < Cft 2 ||VR ||( L 2 (u6h7o)) 3x3. 

The constants depend only on doj. 

Secondly, we set Mh the maximum of the function (i*i,i* 2 ) h-» 1 1 1 13 — R^(i*i, i* 2 )||| into the closed set 
{(1*1,1*2) € oj I ft < ciist((i*i,i* 2 ),7o) < 2ft}, and let (si,s 2 ) be in this closed subset of oj such that 

M h = |||l3-Rfc(«i,s 2 )|||. 
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Applying the Poincare-Wirtinger's inequality in the ball B((si, s 2 ), 4ft) we deduce that 

V(s' 1 ,4) G B((s 1 ,s 2 ),h), IWR'^s'-l, s' 2 ) - R h (s 1 ,s 2 )\\\ < C||VR'||( L 2 (B((si)S2))4h))) 3. 

The constant depends only on the Poincare-Wirtinger's constant for a ball. 
If M h is larger than C||VR'||(i,2( B (( Sl;S2 ) i 4 fe )))3 we have 

nh 2 (M h - C||VR ||(i,2(B(( sllS2 ),4/ l )))3) < IIRft - I 3|||(i,2(B(( Sl , S2 ),/ l )))3 

<ll R ft ~ l3||(Z,2 (u , 3hi7o)) 3X3 < Ck 2 \\\7R || (i 2 (wf .^ 7o)) 3X3 

then, in all the cases we obtain 

M /l <C||VR'|| (L 2 (w6 ^ 7o)) 3x3. 

The constant does not depend on h and R . The above inequalities show that there exists ft such that for 
any h < h 



|R h («i,« 2 )-l3||| <C||VR|| (i 2 ( , 



"6^ )) 3X3 



< 1/8 for any (s\, s 2 ) G oj such that ft < dist((si, s 2 ), 70) < 2ft. 



By definition of R^, that gives |||R/j(si, s 2 ) — I3 j | j < 1/8. 

Finally, for any h < max(h' , ftg) and for any (si, s 2 ) we have 

dist(R h {s u s 2 ),SO(3)) < 1/8. 

Using (C.5) and (C.6) we obtain (recall that || • || 2 is the euclidian norm in 

V(si,s 2 ) € cv, -^(si,s 2 )-K h (si,s 2 )t a (si,s 2 ) 2 < Cft||</>|| (W 2, c » (w)) 3+C(||V||(B2( W 3 ft )) 3 +ll R -'ll(H 1 (w3h)) 3x3 ) 

where oj^h — {{si,s 2 ) e M. 2 \ dist((si,s 2 ), duj) < 3ft}. 
We have 

dVh vt + n 1 \( dV 
^--^t a = (l-^)(- 



^-R^ = (1 -^(^-R^) + ^(,-V;). 



Thanks to the above inequality, (C.4) and again the estimate of |||R/j— 13 1 1 1 in the edge strip h < dist((si, s 2 ), 70) < 
2ft we obtain for all (si, s 2 ) e iv 

^^(si,s 2 ) - R, l (si,s 2 )t Q (si,s 2 ) 

OSa 2 

<C(ft||0||( l y2,oc( w ))3 + ||V||( H ;2 (w3h)) 3 + ||R ||(Bl(u, 3h ))3X3 V||(B2 (u , 5h7o)) 3X3). 

The same argument as above imply that there exists fto < max(/i , ft ) such that for any < ft < ho and for 
any (si, s 2 ) S w we have 



(C.ll) 



dV h 

ds a 



(si, s 2 ) - R/i(si, s 2 )t a (si, s 2 ) 



< 



1 



From (C.4), (C.5), (C.6) and (C.7) there exists a positive constant C which does not depend on h such that 

C 



(C12) 



R 



+ I|R||(h 1 (w)) 3 >< 3 }- 
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Now we can choose h in term of S. We set 



h = 6S, Se (0, So] 

and we fixed 9 in order to have h < ho and to obtain the right hand side in (C.12) less than — _ , ( c x is 

4v2c 1 <!) 

given by (2.3)). It is well-known that there exists a sequence {Vs) s ^ 5o ] such that (Vj, R5, Vj) € B<5, 70 and 
satisfying the convergences in (C.l) and the estimate in (C.3). □ 
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